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Abstract 

We give a deterministic 0(/in^'’'^/^)-time (2/i)-approximation nonadap¬ 
tive algorithm for 1-median selection in n-point metric spaces, where h S 
Z"*" \ {1} is arbitrary. Our proof generalizes that of Chang [2]. 


1 Introduction 


A metric space (M, d) is a nonempty set M endowed with a function d\ MxM —)■ 
[0, cxD) such that for all x, y, z & M, 

• d{x, ?/) = 0 if and only if x = y, 

• d{x,y) = d{y,x), and 

• d{x,y) + d{y,z) > d{x,z) (triangle inequality). 

The METRIC 1-MEDlAN problem asks for a point in an n-point metric space (M, d) 
with the minimum average distance to other points. For c > 1, a point p G M is 
said to be c-approximate for metric 1-median if 


d(p,x) < c ■ mm > d 

peM ^ 

x&M x&M 

An algorithm for METRIC 1-MEDIAN is nonadaptive if the sequence of distances 
that it inspects depends only on M but not on d. Because there are n{n — l)/2 
nonzero distances, “sublinear-time” will mean “o(n^)-time.” 
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Indyk [5, 6] shows that metric 1-median has a Monte-Carlo 0(n/e^)-time 
(1 -|- e)-approximation algorithm for each e > 0. In where D > 1, metric 
1-MEDlAN has a Monte-Carlo 0(2P°h(iA)/))_tinie (l-l-e)-approximation algorithm 
for all e > 0 [7]. Many other algorithms are known for fc-median selection [1, 4, 7]. 

For example, Guha et ah [4] give a deterministic, 0(n^’'''^)-time, 0(n^)-space, 
2'^(^/^)-approximation and one-pass algorithm as well as a Monte-Carlo algorithm 
for fc-median selection in metric spaces, where e > 0. 

We show that metric 1-median has a deterministic 0(/m^+^/^)-time (2h)- 
approximation nonadaptive algorithm for all h G Z’*' \ {1}, generalizing the fol¬ 
lowing theorems: 

Theorem 1 ([2]). Metric 1-median has a deterministic 0{n^'^)-time A-approximation 
nonadaptive algorithm. 

Theorem 2 ([8]). For each h G Z+\{1}, metric 1-median has a deterministic 
0{hn^^^^^)-time {2h)-approximation (adaptive) algorithm.^ 

When n is a perfect square and h = 2, our proof is equivalent to that of 
Theorem 1 [2], Chang [3] shows that metric 1-median has no deterministic 
o(?7,^)-query (4 — f2(l))-approximation algorithms (where an algorithm’s query 
complexity is the number of distances that it inspects). So the approximation 
ratio of 4 in Theorem 1 cannot be improved to a smaller constant. 

2 Our algorithm 

Let ({0, 1,..., n — 1}, d) be a metric space, h G \ {1} and t = . For all 

j G {0,1,..., n — 1}, denote the (unique) f-ary representation of j by 

(Sh-l(j), Sh-2(j), • • •, So(j)) G {0,1,..., f - 1}'", 


i.e.. 


h-i 


seU) ■ = j- 


( 1 ) 


For all i, j G {0,1,..., n — 1}, 


h—1 / h—1 h—1 \ 

d{i^i-\- j mod n) = d\i+ st{j) ' rnod n, f -|- E Si{j) ■ t^ mod n \ . (2) 
k =0 V e=h—k £=h—l—k / 

By convention, empty sums vanish, e.g., = 0. 

Mhe time complexity of 0{hn^'^^F) is originally presented as because h is in¬ 

dependent of n. We include the 0{h) factor, which is implicit in the original proof, for ease of 
comparison. 
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Lemma 3. For all i, j E {0,1,... ,n — 1}, 

d{i,i + j mod n) < d{i,i + j mod n) 
Proof. By equation (2) and the triangle inequality for d. 


h-l 


d{i,i + j mod n) > d \ i,i + S£{j) ■ mod 


n 


i=0 


This and equation (1) complete the proof. 
Lemma 4. For all a G {0,1,..., n — 1} with 


n—1 


n—1 


d{a,a + j mod n) = min d{i,i + j mod n), 

j=0 j=0 


we have 


n—1 


J2d{a,j) < 2h-min ^ d{i,j). 

j=0 j=0 


n—1 


Proof. Let it be a uniformly random element of {0,1,..., n — 1} and 


n-l 


n—1 


It J. ■ ^ 

i' = argmin d {i,j), 


i=0 


j=0 


breaking ties arbitrarily. It is easy to see that 


n—1 


n—1 


J2d{a,j) = J2 d (a, a + j mod n). 

j=0 j=0 


Furthermore, 


n—1 


Lemma 3 


n—1 


d (a, a + j mod n) < d{a,a + j mod n) 

j=0 j=0 

' n—1 

d {u, u + j mod n) 


equation (3) 

< E 


j=0 


By equation (2), 


E 


= E 


n—1 


d (u, u + j mod n) 


□ 


(3) 


(4) 


(5) 


j=0 

n—1 h—1 / h—1 h—1 

EE d j It + se{j) ■ mod n,u + S£{j) ■ mod 

j=0 k=0 V £=h—k i=h—l—k 


n 
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Finally, 

n—1 h—1 


E 


< E 


h-l 


h-l 


EE c? I “t" ^ ^ ^ mod. Ti^u -\- ^ ^ t mod 


n 


j=0 k=0 
n—1 h—1 


i=h-k 

h-l 


l=h—l—k 


h-l 


EE d ii',u+ ^ Si{j) ■ mod n \ + d\i\u+ ^ mod 

j =0 k =0 \ i=h-k J V e=h-l-k 

/ h-l \ 


n 


d\i',u+ Si{j) ■ mod 


n 


l=h-k 


+ E 


h-l 


d\i',u+ si{j) • mod 


t=h—l—k 


n—1 h—1 

= EEe 

j=D k=0 
n—1 h—1 

= Z] Z] +E[d(i',-u)]) 

j =0 fc =0 

= 2nh • E [ d (^^ 'u) ] , 

where the inequality follows from the triangle inequality for d, and the second-to- 
last equality is true because u + ^ distributes uniformly at 

random over {0,1,..., n — 1} for any j G {0,1,..., n — 1} and A; G {0,1,, h}. 
Inequalities (5)~(6) imply inequality (4). □ 

For a predicate P, let x[P] = 1 if P is true and x[P] =0 otherwise. Dehne 
s'^_ 2 i • • •) Sg) G {0,1,..., t — 1}^ to be the t-ary representation of n — 1. So 
J2r=o = n — 1. For i G {0,1,..., n — 1} and m G {0,1,..., h — 1}, 


n 


z ^ 


— ! ,...,So— 0 L?’ — 0 




r=0 


d\i+ Si ■ mod n, i si ■ mod 


n 


(7) 


k=0 


i=m+l—k 


l=m—k 


t-1 


g{i,m) = E E d I i + S£ ■ mod n,i + S£ ■ mod n j . (8) 


—1 —0 k —0 


£=yrL-\-l — k 


i=m—k 


Clearly, 


/(bO) = E d{i,i + So mod n ), 
* 0=0 
t-i 

d(b0) = E d{i,i + So mod n). 


(9) 

( 10 ) 


* 0=0 


Lemma 5. For all i & {0,1,... ,n — 1}, 


n—1 


f {i,h-l) = ^ d{i,i+j mod 

j=0 


n) 
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Proof. As Ylr=o s'j. ■ P = n - 1, 


t-i 


X 


h-1 


Sr ■ P < n — 1 

r=0 

Yd[^+ Y Si ■ mod n, i + 'Y^ Si ■ mod n j . (11) 


h-l 


h-1 


h-l 


k=0 


£=h-k 


£=h—l—k 


By the existence and uniqueness of a t-ary representation of each j G {0,1,..., n— 

1 }, 

n—1 h—1 / h—1 h—1 \ 

EE c? I z “h ^ ^ ^ mod. Ti^ % “h ^ ^ ^ mod ti | 

j=0 k=0 \ l=h—k i=h—l—k / 

h-l 

'Y, Sr ■ P < n — 1 


Y ^ 

Sh-l,Sh-2,-"iSO=0 


r=0 


h-1 


h-1 


h-1 


■ E‘*r+ E Si ■ mod n,i + si ■ mod n 

k=D \ i=h—k i=h—l—k / 

Equations (2) and (11)-(12) complete the proof. 

Lemma 6. For all i & {0,1,... ,n — 1} and m G (1, 2,... ,h — 1}, 

t-i 

g{i,m) = " Y ^ d (z, i + Sm ■ mod n) 


( 12 ) 

□ 


—0 

t-1 


"Y^ 9 {'l P Sm ■ mod n,m — 1). 


d{i,i + Sm ■ P^ mod n) 


Sm —0 

Proof. By equation (8), 
g{i,m) 

t-i t-i / 

= E E 

Sm—O Sm-l,Sm-2',---,SO—0 

m— 1 / 

+ ^ d i i + Sm ■ ^ S£ ■ P mod n,i + Sm • ^ se ■ P mod ^ ) ) j 

A;—0 \ i—m—k i—m—l — k / / 

g [i + Sm ■ mod n, TO — 1) 

t—1 m—1 / 771—1 771 — 1 

= ^ ^ d\i + Sm-P^+ n,i + Sm ■ P^ + ^ sg ■ P mod ; 


771— 1 


771— 1 


Srn - l , Sm - 2 ,---, SO =0 fc =0 


i—m — k 


i—m—l — k 


for Sm G (0,1,..., t — 1}. Furthermore, 


t-i 


t-i 


t-i 


El El d {i,i + Sm ■ mod n) = El d {i,i + Sm ■ mod n). 

Sm, _ 1,5771— 2,. 


Sm,— 0 
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Lemma 7. For all i E {0,1,... ,n — 1} and m E {1,2,... ,h — 1}, 


□ 


( m-l \ 

1 + f \ d{i,i + s'^ ■ mod n) 

+ ^ d{i,i + Sm ■ mod n) 

-5 m— 0 

+ / (i + mod n,m — l) 

+ 5f (i + Sm ■ mod n, m — 1). 

Sm— 0 

Proof. Observe the following for all Sm, Sm- 1, ..So G {0,1,..., t — 1}: 

(i) If Sm = s'^, then Y.T=o Sr ■ F < Y.T=o ^'r ■ if and only if Y.7=o Sr ■ F < 

y^m-l / y.r. 

Z^r=0 ^ 5 

(ii) If < s'^, then ^ Sr ■ F < X;™ o 4 ' f"'; 

(hi) If s^ > s'^, then X)™ o Sr ■ F > X)™ o 4 ' f"'- 
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We have 


equation (7) 


item (iii) 


item (i) 


item (ii) 


d{i,i + Sm ■ W mod n) 


t-1 t-1 

E E X 

S^_i,S^_2,...,So=0 

m—1 / 

+ d j i + Sm • W + St ■ mod n, i + • W + st ■ mod ' 


(13) 




_r—0 
m— 1 


r=0 


fe=0 


£—m—k 


i—m — l — k 


t-1 

E E 

5^71=0 Sm-l,Srn,_2,..-,So=0 




+x 


(Sm = sL) X ( ^ Sr • i’' < XI 

r—0 / 

d{i,i + Sm ■ W mod n) 


m— 1 





(Sm < S(„) A I X Sr • i’' < X 

\r—0 r—0 

1 / m— 1 

+ X d I * + Sm • W + X n, z + Sm • W + X s^ • mod n 

i—m—k ^—m — l — l 

/m—1 m —1 

(Sm = S'm) A I X Sr • i’' < X '®r ’ 


E E 

5^71=0 S77i_i,S7„,_2,...,So=0 


m— 1 


, r—0 


r—0 


+x 


fe=0 


l—m—k 



d{i,i Sm ' mod n) 


t-1 


(Sm < S'm) A ( X Sr • i’' < X ^r ’ 

\r—0 r—0 

1 / m— 1 

+ X d I z + Sm ■ W + X Sf • mod n, z + Sm • W + X s^ • mod ' 

i—k i.—m — 1 — i 

( m—1 m —1 

X Sr • i’' < X '®r • i’' 

r—0 r—0 

+X [sm < Sm] j • l^d (z, z + Sm • w mod n) 

m—1 / m—1 m—1 

+ X d I z + Sm ■ W + X s^ • mod rz, z + Sm • W + X s^ • mod ' 


E E U 

Sm—^ S77i_i,S7„_2,...,So=0 \ 


fe=0 


£—m—k 


£—m — l — k 
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By equation (7), 

/ (i + • t™' mod n, m — l) 


t-i 

E 


X 


Sm — 1 ?S 77 j ,_2 — 0 L T — 0 


m—1 


m—1 


Y, Sr-f <Y, ^r- t" 


r=0 


m—1 


m—1 


m—1 


£=m—k 


d i i + s'jn ■ + Y^ ri,i + s'm ■ + Y^ ^ 

k=0 V 

t-1 

Sm —0 Sm — 1 ,Srn, —2 — 0 


^m—1 


i=m—l—k 
m—1 ' 


(Sm = s'm) A ^ Sr • t"” < J]] s'r ' f 


. r=0 


r=0 


m—1 


m—1 


m—1 


Y, d[i + Sm-f^+ Y^ n, i + Sm • + Yj 

/c^O V i=m—k i=m—l—k 

By equation (8), 

sin-1 

Y^ 9 {i + Sm ■ mod n,m — 1) 

Sm—0 

Sm t-1 

= X] X] x[sm<sY 

Sm=0 Sm — 1 5^m, —2 )”->'50~0 


n 


m—1 


m—1 


m—1 


Y, d[i + Sm-f^+ Y^ n, i + Sm • ^ S£ • mod 


n 


A:=0 


i=m—k 


i=m—l—k 


Because each number in {0,1,..., YY=o ^'r ' t)e written uniquely as 

Y1T=0^ ■Sr ■ t\ where Sm-l, Sm-2, ..., So e {0,1,..., t - 1}, 


i-l 


[ ^m, — S 


m—1 m—1 

r=0 


, r=0 


Sm — 0 Sm — l:^m — 2f‘-’^^0 —^ 

d (i, i + Sm ■ mod n) 

sjn / m-1 \ 

Y h- + X] '®r ■ ■ xlsm = s'^] ■ d(i,i + Sm -t"" mod n) 


Sm —0 


r=0 


m—1 


1 + Yj Sr I d{i,i + mod n). 


r=0 


Finally, 


t-i 


E E x[sm < s'^]- d{i,i + Sm-d^ mod n) 

-5m—0 Sm — 1 j-^m —2 V - 0 

^m-1 

'Y^ d {i,i + Sm ■ t™' mod n). 


Sm —0 


( 14 ) 








1: t <(— ; 

2: Find the t-ary representation of n — 1, denoted s'f^_ 2 , ■ ■ ■, Sq) G 

3: for i G {0,1,..., n — 1} do 
4: /[^][0] ^ E^o=o + modn); 

5: ^[i][0] ^ E*.o=o * + So mod n); 

6: end for 

7; for m = 1 up to h — 1 do 
8: for i G {0,1,..., n — 1} do 

9: f\i][m] ^ (1 + EE"o^ K ■ f) d{h i + mod n); 

10: /[*][m] ^ /[*][m] + t”* Es^=o * + Sm ■ mod n); 

11: /[*][m] ^ /[*][m] + f[i + s'^ -mod n][m — 1]; 

12: f[i]['m] ^ f[i][m\ + g[i + Sm ■ mod n][m - 1]; 

13: g[i] [m] •(— E1~Lo ^ ■ t"^ mod n); 

14: g[{\ [m] ^ g[{\ [m] + E1~Lo 9[i + Sm ■ t"' mod n] [m - 1]; 

15: end for 

16: end for 

17: Output argmiuEo^ /[*][^ ~ I]; breaking ties arbitrarily; 

Figure 1: Algorithm find-median with input a metric space ({0,1,..., n — 1}, d) 
and h G Z’*' \ {!}. 

Equations (13)-(14) complete the proof. □ 

Lemma 8. Algorithm find-median in Fig. 1 is {2h)-approximate for metric 1- 
MEDIAN. 

Proof. By equations (9)-(10), lines 4-5 of find-median compute f{i, 0) and g{i, 0). 
Then, by Lemmas 6-7, f{i,m) and g{i, m) can be found by dynamic programming 
as in lines 9-14. So line 17 outputs argminEo^ /(i, h — 1), which equals 


It ± ^^ ~ 

argmin y d (z, i -|- j mod n) 
i=o 

j=0 

by Lemma 5. Now Lemma 4 gives the approximation ratio of 2h. □ 

We now state our main theorem. 

Theorem 9. Metric 1-median has a deterministic 0{hn^^^/^)-time {2h)-approximation 
nonadaptive algorithm for each h G \ {!}. 

Proof. Clearly, find-median is deterministic and nonadaptive. Furthermore, it is 
(2h)-approximate for metric 1-median by Lemma 8. As s' < f — 1 for all 
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i E {0,1,..., /i — 1}, the loop in lines 3-6 of find-median takes 0{nt) time. By 
precomputing f and X]r=o all i G {0,1,..., h — 1}, each iteration of the 

loop in lines 8-15 takes 0{t) time. □ 
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